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ABSTRACT 

We study the collision rates and velocities for point-particles of different sizes in turbulent flows. We 
construct fits for the collision rates at specified velocities (effectively a collisional velocity probability 
distribution) for particle stopping time ratios up to four; already by that point the collisional partners 
are very poorly correlated and so the results should be robust for even larger stopping time ratios. 
Significantly, we find that while particles of very different masses have approximately Maxwellian 
collisional statistics, as the mass ratio shrinks the distribution changes dramatically. At small stopping 
time ratios, the collisional partners are highly correlated and we find a population of high number 
density (clustered), low relative- velocity particle pairs. Unlike in the case of identical stopping time 
collisional partners, this low relative- velocity clustered population is collisional, but the clustering is 
barely adequate to trigger bulk effects such as the streaming instability. We conclude our analysis 
by constructing a master fit to the collisional statistics as a function only of the stopping time ratio. 
Together with our previous work for identical stopping time particle pairs, this provides a recipe for 
including collisional velocity probability distributions in dust coagulation models for protoplanetary 
disks. We also include our recipe for determining particle collisional diagnostics from numerical 
simulations. 

Subject headings: turbulence - planets and satellites: formation 



1. INTRODUCTION 

Turbulent collisions between small dust grains is be- 
lieved to lead to sticking and growth of the solids in 
protoplanetary disks, a crucial early step in planet for- 
mation. However, as dust grains grow in size, they 
respond more slowly to the gas motions. This allows 
two neighboring dust grains to feel the effect of larger 
scale, higher velocity gas motions differently, leading 
to higher collisional speeds. Once the particles collide 
with large velocities, these collisions result in growth- 
neutral bouncing events and, when even faster, frag- 
mentation. While the precise speed of turbulent mo- 
tions in protoplanetary disks and the material strength 
of the dust is unknown, the requirement that turbu- 
lence be adequately strong to explain observed accre tion 
rates (Shakura fc Sunvaevlll973t iRussell et al.ll2006j ) im- 
poses peak turbulent velocities noticeably above lOm/s, 
which itself is well a bove the dust fragmentation thresh- 
old of perhaps 1 m/s (|Blum fc Wurmll2008t [Guttler et al.l 
2010). It follows that there is a maximum grain size be- 
yond which equal-mass collisions will usually be too fast 
for the participants to survive; and so such collisions lead 
to fragmentation rather than further dust grai n growth 
(jDullemond fc Dominikl [20051 : iZsom et al.ll2010t ). 

Historically, most of the research into turbulence in- 
duced dust-dust collisions in astrophysics was spent an- 
alytically determining a single characteristic collisional 
velocity scale for two dust grains of specifie d frictional 
stopp i ng times with re s pect t o the gas dVolk et al.l 
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Given the large turbulent velocities required to achieve 
significant and observed accretion flows in protoplane- 
tary disks, this has led to the conclusion that bounc- 
ing and fragmentation barriers to collisional dust grain 
growth are inevitable. However, some recent work has 
instead focussed on the effect of including not merely a 
single collisional velocity scale, but rath er an entire col- 
lision al velocity probability distribution (Okuzum i et al.l 
1201 If ). Imposing Maxwellian distributions around the 
previously estimated characteristic collisional velocities, 
they found that small numbers of lucky grains who 
have multiple consecutive low-velocity encounters could 
grow significantly past the bouncing and even frag- 
ment ation barriers ( Windm ark et al1l2012t iGaraud et al.l 
2012). This growth naturally peters out as large lucky 
grains are rare, and eventually unlucky collisions de- 
stroy them. However, if these lucky grains become large 
enough compared to the bulk population, which remains 
trapped by the bouncing barrier or ground up by the 
fragmentation barrier, then they can grow even as the re- 
sult of high velocity collisions with much smaller grains in 
a process is referred to as mass transfer: they are large 
enough to survive the hits. The grain growth process 
then is dominated by the sweep up of small grains by a 
lucky few large dust grains, with equal mass collisions be- 
tween the large grains being destructive but, fortunately, 
vanishingly rare due to their small nu mber density. 

In a previous paper (|Hubbardll2012l . henceforth Paper 
I), we looked at the collisional velocity distribution for 
a monodisperse population of dust grains (i.e. a popula- 
tion with the same gas-dust drag parameters) , and found 
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a distribution that was both lower in velocity than previ- 
ously expected by a factor of approximately five, and very 
non-Maxwellian. Indeed, in that case the collisional ve- 
locity distribution had a pure exponentially decaying tail, 
and the particle pair velocity probability distribution ap- 
peared to diverge at low relative velocity. In this paper 
we extend that analysis to include polydisperse popula- 
tions (collisions between grains of differing sizes), going 
from gas-dust frictional stopping time ratios of 17/16 to 
4. The stopping time scales with the grain radius and so 
the cube root of the mass, so our largest stopping time 
ratio is equivalent to a mass ratio of 64, large enough 
for mass transfer to occur. It is also large enough that 
we clearly approach an asymptotic limit as the two dust 
grains become uncorrelated. 

Our goals are two-fold. Firstly, we simply wish to un- 
derstand any features and behavior of the collision veloc- 
ity probability distribution as there is no analytical treat- 
ment that can handle long time correlations of similar 
sized dust grains, and such correlations were a dominant 
feature in the monodisperse case. Secondly, we wish to 
develop a master fit for the collision velocity probability 
distribution as a function only of the relative size of the 
dust grains. Such a fit is of clear utility for any dust grain 
collisi on gr owth model su c h as u sed in IWindmar k et ahl 
(12012ft and iGaraud et all (|2012D . We should note that 
our master fit only applies in the limit of the dust grains 
being well contained within the turbulent cascade, as will 
be discussed below. 

In Section 2, we discuss our turbulence model and how 
we extract the collisional diagnostics, going into greater 
detail on the latter in the appendix. In Section 3, we 
perform initial analysis of our results, extracting cluster- 
ing information as well as the most common turbulent 
collisional speed. In Section 4 we examine the shapes of 
the collisional velocity probability distributions and cre- 
ate fits to them which can be used in further modeling 
of dust coagulation. We conclude in Section 5. In the 
Appendix we include a detailed procedure for extracting 
dust collisional diagnostics from numerical simulations. 

2. METHODS 

As in Paper I, we evolve particles in a semi-analytic 
model for turbulence to extract collisional velocity prob- 
ability distributions. We do this for a range of stopping- 
time ratios, constructing fits to the collisional velocity 
probability distribution for each ratio, with a greater em- 
phasis on identifying features in the distributions than 
pure quality of fit. Finally, we construct a fit to the fits 
to arrive at a final fit formula of the collisional velocity 
distribution with a small set of tunable parameters that 
can eas ily to used in further pa rticle coagulation studies 
such as IWindmark et ahl (|2012h . 

2.1. Particles 

We insert 10 6 particles into the synthetic turbulence 
field. Writing u for a dust particle's velocity and v for 
the gas velocity at the particle's position, the particle's 
equation of motion is 

du u — v 

~dt = ~ ( > 

where r is the particle's stopping time. Note that r is 
proportional to dust grain radius for constant grain mate- 



rial density. All our simulations split the particles evenly 
between two species that differ only through their stop- 
ping times. We refer to the larger stopping time as t\ 
and the smaller as T2- Further, we define 

e = H - l > (2) 

as the measure of the stopping time ratio. 

The collisional properties of a dust grain pair in in- 
compressible turbulence are expected to scale with the 
stopping time of the largest grain because the shorter 
stopping time grains are more tightly bound to the gas 
and particle collisions require deviations from the incom- 
pressible gas velocity. Further, dust is expected to couple 
(for collisional purposes) most strongly with turbulence 
with turnover times it = r. This is because if the turbu- 
lence has a turnover time t t S> r the dust grains will be 
completely entrained by the flow on a time scale t and 
will barely deviate from the motion of the gas; while if 
the turbulence turnover time it <C t the dust grains will 
not have time to respond to the gas motion. Turbulent 
motions with timescales it = r are therefore Goldilocks 
motions, both varying rapidly enough that the particle 
motion can deviate from that of the gas, and varying 
slowly enough that the particle motion is measurably af- 
fected. 

We therefore denote u p and k p as the velocity scale 
and wavenumber of the turbulence with turnover time 
l/kpU p = t\ = t t (kp). We denote the associated veloc- 
ity scale as u p (a dust velocity) rather than v p (a gas 
velocity) because this is the velocity scale we will use 
to normalize our dust collisional velocities. When k p is 
not actually included in our turbulence model (which is 
discrete in k), u p is determined by interpolating the Kol- 
mogorov spectrum. 

We consider modest e < 3 due to both practical con- 
siderations (numerical resources) and because the larger 
e case is expected to be straightforward. That is be- 
cause once the stopping time ratio becomes large, the 
two populations of dust grains become uncorrelated, and, 
from the perspective of the particles with stopping time 
T\ T2 , the grains with stopping time t-i are completely 
bound to the gas motion and so indistinguishable from 
the gas. This means that the collisional behavior will 
have a well defined asymptotic behavior, which we find. 

2.2. Synthetic turbulence - Motivation 

Since iKolmogorovl (|1941h , turbulence has been under- 
stood as a cascade of kinetic energy from one lengths 
scale to a shorter one. The largest length scale is the 
one associated with the driving source of the turbulence, 



TABLE 1 

Runs 



Run 


€ 


Tl 


T2 


tls/ri 


T2/tss 


A 


0.0625 


2.01875 


1.9 


5.0 


7.7 


B 


0.125 


1.96875 


1.75 


5.1 


7.1 


C 


0.25 


2 


1.6 


5 


6.5 


D 


0.375 


2.3375 


1.7 


4.3 


6.9 


E 


0.5 


2.25 


1.5 


4.4 


6.0 


F 


0.9 


2.47 


1.3 


4.0 


5.2 


G 


1.25 


2.7 


1.2 


3.7 


1.8 


H 


1.4 


2.88 


1.2 


3.5 


1.8 


1 


2 


3 


1 


3.3 


4.0 


J 


2.5 


3.5 


1 


2.9 


4.0 


K 


3 


4 


1 


2.5 


4.0 
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while the smallest length scale is determined by the vis- 
cous dissipation of the turbulent energy into heat. Be- 
tween those two loosely defined scales lies the inertial 
range. In the inertial range the turbulence is scale free, 
independent of the details of the forcing and of the dis- 
sipation: the only parameter is the energy cascade rate. 
The Kolmogorov power spectrum of v(k) oc A: -1 / 3 de- 
rived from the simplest dimensional analysis of the prob- 
lem fits experimental data well. Significantly, in proto- 
planetary disks, the inertial range is believed to be large 
in that the ratio of the time scale ti s associated with the 
largest scale motions to the time scale t ss associated with 
the smallest scale motions is large: ti s ^> t ss . 

From the hypothesis that particle-particle collisional 
properties arc determined by turbulence with t t ~ t\, 
and the consideration of modest e it follows that there 
are five regimes for turbulence induced particle-particle 
collisions, which can be describe by relating r to the 
bounding turbulent time scales, U s and t ss . The regimes 
are: 

• Particle stopping times much longer than the 
longest turbulent timescale, or r 3> U s : This 
regime may be accessible to direct simulation of the 
turbulence using Navier-Stokes with current nu- 
merical resources, although the particles will move 
significant distances so the numerical domains will 
have to be large. 

• Particle stopping times comparable to the longest 
turbulent timescale, orr~ ti s : This regime should 
be accessible to direct numerical simulation, al- 
though the applicability to protoplanetary disks 
may be complicated by rapid particle drift through 
the disk, and the results may depend sensitively on 
the turbulent driving mechanism. 

• Particle stopping times comparable to the shortest 
turbulent time scale, or t ~ t ss : This regime is 
numerically accessible for hydrodynamical simula- 
tions, although the dissipation must be well cap- 
tured. 

• Particle stopping times much shorter than the 
shortest turbulent time scale r <C i ss : This regime 
is numerically accessible. 

• Particle stopping times much longer than the short- 
est turbulent time scale, but also much shorter than 
the longest turbulent time scale, or ti s 3> r 3> t ss : 
This regime is where the bouncing/fragmentation 
barriers are believed to lie, but is numerically inac- 
cessible because such an inertial regime is beyond 
current numerical resources. 

We are interested in the last regime, because full col- 
lisional velocity probability distributions are interesting 
precisely when collisions at different speeds have differ- 
ent outcomes. This regime is, however, not one that can 
currently be achieved through direct numerical simula- 
tion. This is our motivation for using semi-analytical 
synthetic turbulence: the synthetic turbulent cascade 
perfectly matches the Kolmogorov spectrum that fits ex- 
periment well while having a large inertial range and be- 
ing numerically feasible. 



2.3. Synthetic turbulence - Model 

Our turbulence field uses the same approach as in Pa- 
per I: a synthetic turbulence field that matches the power 
spectrum of the ine rtial range of Kolmogorov turbulence 
(|Kolmogorov| [r941') . In other words, while the absolute 
scale of the kinetic energy is arbitrary, the relative en- 
ergy in the motions at two given scales exactly matches 
that of perfect Kolmogorov turbulence, which matches 
experiment well. 

In this model, the turbulence is binned logarithmically 
in wave- vector space. A region of wave- vector space with 
|fe| € fc re f ± 5k is a spherical shell, so we refer to the bins 
as such, and the shell m is associated with |fc| = 2 m . 
Each shell m is assigned a trio of wave-vectors k mn of 
the appropriate magnitude. The gas velocity associated 
with each wave- vector k mn is along a velocity unit vector 
v mn _L k mn , so that the final velocity field is incompress- 
ible. 

In Paper I we used multiple methods for introducing 
time variation in the velocity field, including rotating the 
projection of the energy in shell m onto the trio of wave- 
vectors and varying the phase factor. In this paper, we 
will use the rotation of the projection, along with a slow 
variation in phase. The rotation of the projection occurs 
as a random walk of length 2n over an eddy turnover 
time t m — l/v m k m , and allows the velocity at a specific 
position to rotate in space, avoiding steady velocity pat- 
terns. The phase variation takes the form of a random 
walk of the phase m „ of length 7r over 3 eddy turnover 
times. The final gas velocity field is 

m,n— 8,3 

V(x, t) = ^2 V2a mn (t)v m V mn COS [k m n ■ X + <p mn (t)] , 
m,n— 0,1 

(3) 

where v m is the velocity associated with the shell m and 
a m (t) is a unit vector that rotates as a random walk 
over the turbulent time scale t m . We set vo = 0.1 in 
code units. We refer the reader to Paper I for a more 
complete description and derivation. We will discuss the 
need for the phase variation in Section [2.51 

In Paper I, we found that the relative particle/gas mo- 
tion in our model had reasonable statistics compared 
to predictions (note that particle/gas statistics are not 
yet well constrained as observed clustering i mmediately 
impli es non-trivial particle/gas correlations, IPan et all 
[201 If ) ; we also found that the particles were transported 
away from their positions according to a random walk, 
appropriate for turbulent particle diffusion. In this work, 
we also find in Section l3.1l that the relative motion of par- 
ticles with large stopping time ratios do approach naive 
expectations of uncorrelated motion. 

In all simulations in this paper, we include shells with 
m ranging from to 8 (k m = 2 m extending from 1 
through 256). The smallest wavenumber, k m= o = 1, cor- 
responds to our box scale L — 2ir. We set ti s = 10 and 
t ss = 0.25 in code units, where here Is and ss refer to the 
largest and smallest eddy scales included in the simula- 
tion, or k = 1, 256. Extending the Kolmogorov spectrum 
then, there is a k value associated with any ti s > t t > t ss 
such that t t = k~ 1 v(k)~ 1 , although that k value gener- 
ally does not fall on one of the values associated with our 
discrete shells. In those cases we simply extrapolate the 
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Kolmogorov spectrum. The large value of ti s /t ss — 40 
allows us to fit a particle pair with stopping time ratio of 
t\/t2 — 4 while maintaining at least U s /t\ ~ Tijt ss ~ 3. 
Our particle pairs are then contained within the effective 
turbulent cascade. We list our stopping time choices in 
Table [U 

2.4. Dust snap-shots 

Our particles are initialized with a random position 
and zero velocity. In Paper I, we found that it took a 
significant time for the particle collisional diagnostics to 
achieve a steady state, in part due to the formation of 
dense clusters of highly correlated particles. Following 
that result, we take snap-shots every turbulent turnover 
time for the largest eddy, starting at the 20th turnover. 
We then average over the snap-shots, identifying small 
separation particle pairs and binning them in relative 
velocity (bin width 6u). 

For a detailed description of our procedure for extract- 
ing dust collisional diagnostics from this data, see the 
appendix. In brief, we use the binned data to deter- 
mine N(R,u, j), which is the number of particle pairs 
(a, b) with separation \x a — Xf,\ < R, and relative ve- 
locity \u a — Ub\ G u ± Su/2. The parameter j repre- 
sents whether the particles have r = Tb = t\ (j = 1), 
T a = n = r 2 (j = 2) or r a 7^ T& (j = 3). Note that the 
largest R we consider is Ro with Rok m= o — 0.02. Com- 
paring with the minimum velocity limit u m i n from the 
appendix, we note that the minimum relative velocity we 
can consider is 



0.2- 



Rn 



I p lip 



0.2^ > 7 x 10- 3 . (4) 
L Lp 

Our analysis remains are safely above this floor. 

From N we derive the fractional density enhancement 
over background 



p(R,u,j) = N(R,u,j)/ 



n 47ri? J 



2 3 L 3 



/Su, (5) 



where n is the number of particles, L? is the box vol- 
ume and we treat p as a smooth function. For both ./V 
and p, a dropped j implies that we are considering only 
collisions between particles of differing stopping times 
(j = 3) while a dropped u implies summation over all ve- 
locity bins (N) or integration over velocity (p). A value 
of p{R,j) = 1 is consistent with randomly distributed 
particles while a value of p(R,j) > 1 implies clustering 
on lengthscales of R or less. 

As discussed in Paper I, astrophysical particles are gen- 
erally point particles in that their size is much smaller 
than that of relevant gas dynamical length scales. Ac- 
cordingly, we will attempt to determine the limit 



P(u,j) = Urn p{R,u,j). 



(6) 



The importance of this step can be understood through 
a consideration of the "cold" population found in Paper 
I. This population was large (it dominated the particle 
pair counts) but interacted with itself only with velocity 
w C oid oc R. The net result was that the cold popula- 
tion interacted at such low velocities that it provided 
no actual collisions despite representing the bulk of the 
particle pairs. 




Fig. 1. — Values of p(u,j) from Run K with the time vary- 
ing phase turned on (red/dashed, yellow/dash-dotted) and off 
(black/solid, blue/dotted). The red/dashed and black/solid curves 
are for R = Ro, j ' = 1 while the green/dash-dotted and blue/dotted 
curves are for R = Rq/2, j = 2. Without the phase variation, we 
see an unphysical high u tail which is not self-similar (different 
dependencies on the appropriately scaled u p ). With the phase 
variation, this behavior is eliminated. 

As long as both particles are well trapped within the 
inertial range of the turbulence, then the problem be- 
comes scale free because the particles are too small to 
interact with the driving scale of the turbulence and too 
large to interact with the dissipative scale. It follows that 
all dust collisional parameters, once scaled to u p , must 
be constructible from e, and knowledge of the power law 
of the turbulent cascade. This means that p(u,j), for 
our fixed turbulent cascade, should be independent of 
the precise values of Ti 2 so long as they are both smaller 
than ti s and t ss and u is normalized to u p . 

2.5. Scale free behavior 

As we have emphasized above, a crucial aspect of being 
well contained within the scale free inertial range of the 
turbulence is that the particle-particle collision velocity 
probability distribution must be scale free. Here we use 
the j — 1,2 data to verify that our method produces 
scale free results. 

Turbulent velocity fields with correlation length I are 
expected to be uncorrelated on lengthscales I ^> i. How- 
ever, Equation ([3]) has infinite correlation lengths, which 
might be expected to set up undesirable correlations be- 
tween particles, which we do in fact find. In Figure [IJ 
which uses data from both Run K and an equivalent 
run with the time- variation of the phase turned off, we 
show p(Rq,u) for j — 1 (black/solid and red/dashed) 
and p(R /2,u) for j = 2 (blue/dotted and green/dash- 
dotted). The black/solid and blue/dotted curves do not 
include the time-varying phase. The x-axis is plotted in 
units of Up for the particle pair, which is the gas velocity 
associated with turbulence with turnover time equal to 
the particle stopping time, different for j = 1,2. This 
means that the velocity scaling factor for any given par- 
ticle pair is the one appropriate to that pair. Run K was 
chosen for illustrative purposes as it uses the largest and 
smallest values for ri, T2- 

The size of Rq , chosen to guarantee a reasonable num- 
ber of pair counts at large relative velocity, is too large 
to derive any actual diagnostic values from the curves 
plotted, but a clear difficulty emerges. When the phase 
variation is turned off, a strong high velocity tail emerges, 
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which is not scale free (i.e. does not scale with u p , the 
turbulent velocity scale that the particles interact with) 
as is shown by the distinct slopes for the two particle pop- 
ulations j = 1, 2. In practice, the tail appears to be con- 
trolled by the largest scale turbulence included in the sys- 
tem. We attribute this problem to the infinite correlation 
length scale of Equation © which allows particle motion 
to be correlated over long distances if the stopping time 
ratio is close enough to 1, while physical turbulence has 
well defined and finite correlation lengths. When we con- 
sider polydisperse cases without time- varying phase, the 
effect remains noticeable for e < 0.25, but larger stopping 
time ratios break the correlations. 

We break this correlation for all e by adding the slow 
phase variation in Equation ([3]). While mutually ap- 
proaching particles do sample the same infinite corre- 
lation length gas velocity field, the nodes (v mn (x) = 0) 
that the two particles see will be in different locations as 
they cross positions of the same k ■ x p value at different 
time, with different phases <f> m n(t)- in Figure [1] we see 
that with the phase variation turned on, the particle pair 
distribution reverts to being scale free as required by the 
turbulent cascade. 

This behavior was not seen in Paper I even though 
there were setups where it should have occurred. This 
is due to Paper Fs low velocity bin cutoff and focus on 
small separations. At a given relative velocity u = au p , 
a particle pair can only close a distance of order T\u — 
a/kp. At low relative velocities (it < u p , the case in Paper 
I), particles can only close separations smaller than the 
appropriate turbulent length scales, so the turbulent field 
we used was appropriate. 

3. EXISTENCE OF LIMITS R -> 

Dust grains in protoplanetary disks, generally less than 
centimeters in size, are much smaller than the smallest 
scale turbulent motions expected, which are larger than 
kilometer scales. Accordingly, actual collisions between 
such dust grains cannot be resolved in a numerical sim- 
ulation. As explained in Section 12.41 we therefore take 
the limit R — > for our diagnostic quantities. For finite 
e (unlike the e — case!), this limit shows well defined 
and finite behavior for both the most common relative 
velocity and the total number of particle pairs. 

3.1. Peak collisional velocity 

In Paper I, we found that the location um of the peak 
of N(R, u, j) in velocity space, for particles with the same 
stopping time, is linear in R, with a intercept of 0. This 
implies the bulk of the particle population has a velocity 
distribution that is differentiable in space (i.e. collision- 
less). Allowing for different stopping times produces a 
different result, as shown in Figure [21 where the location 
of the peak remains finite even in the limit R — > 0. For 
stopping time ratios close to e ~ 0, the location of the 
peak (at R — > 0) is approximately linear in e, while it 
levels out for large ri/ra, as seen in the lower panel of 
Figure [21 This is expected as in that limit, much of the 
relative velocity comes from the relative velocity between 
the larger dust grains and the gas (to which the smaller 
grains are tightly bound). 

The fit in the bottom panel of Figure [21 implies a value 
of um — 1.3u p in the limit e — > oo. While we have 
not yet considered the full probability distribution of the 
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Fig. 2. — Top Panel: um{R), the velocity at which p(R,u) has 
its maximum value when taken as a function of u for a range of 
R (black/solid except for H and I which are black/dashed), along 
with their linear fits (red/dotted). There are finite, well defined 
limits of um as R — ¥ 0. From bottom to top are runs A-K in order 
(H,I are dashed for clarity). The inner cut-off is where the pair 
counts become unreliable for the fitting purposes. Bottom Panel: 
Location of the peak (limit R — y 0) as a function of e. 

collisional velocities, this is of the same order as found 
analytically for this case; which is the limit for which 
theory should be most applicable as th e long-time cor- 
relatio ns are the weakest. In particular, lOrmel fe Cuzzil 
(2007) predict that the RMS-averaged collisional veloc- 
ity should reach y/3u p in this limit, which is comparable 
to our limit for um- We can also expect from Figure [2] 
that collisions will be significantly slower than this value 
for e < 1. As particle mass scales with r 3 , this corre- 
sponds to an already noteworthy mass ratio ~ 10. One 
of the greatest difficulties in growing particles collision- 
ally is the fragmentation barrier, which is mitigated if the 
larger grain is larger enough to survive bombardment. 
While the collisional velocity we find is lower than that 
previously anticipated for e < 1, it is still high enough 
to lead to fragmentation for the largest scale turbulence 
expected in protoplanetary disks. Easing the collision 
velocity between near-equal mass grains should however 
lead to a larger production of "lucky" grains. 

3.2. Clustering 

Turbulence has long been known to expel heavy 
(higher density than the fluid) particles from high vor- 
ticity regions through centrifugal forces and concen- 
trate them in th e high strain regions between vor- 
tices (IMaxevlll987HFessler et al.lll994HCuzzi et alj|2001t 
iToschi fc Bodenschatzl 120091 : iPan et al.l 120111 ). Indeed, 
particles with identical stopping time are strongly clus- 
tered in that the dust grain density perceived by a 
test grain diverges to infinity as that grain considers 
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Fig. 3. — Top Panel: Effective density enhancement p(R) for a 
range of R (black/solid), along with their linear fits (red/dashed). 
As in Figure [2] there are well defined limits as R — > 0, which is not 
the case for e = 0. From top to bottom are runs A-K, in order. 
The inner cut-off is where the pair counts become unreliable for 
the fitting purposes. Bottom Panel: Limit R — > of p(R) as a 
function of e. The divergence is weaker than e — t1,77 . 

ever smaller spheres around itself (|Pan et al.ll201lL Pa- 
per I) . If such clustering persists for large ranges in the 
stopping time ratio, the streaming i nstability or gravi- 
tatio n al collapse could be trigger e d (iGoldreich fc Ward! 
[1971 IJohansen et all 120061 [20071: IShi fc Chiang! 120131) . 
We should note here that the vortices we consider 
are too small to feel any Coriolis forces that exist in 
a protoplanetary accretion disk. Large anti-cyclonic 
vortices in accretion disks are kn own to concentrate 
particles through orbital in teractions) Bar ge fc Sommerial 
ll995HJohansen et al.ll200l , but such effects require both 
rotation (which we do not include) and large enough 
length scales for orbital dynamics to play a major role 
(which we do not consider). 

In Figure [3]we show the clustering, i.e. p(R), as a func- 
tion of R and e. In the case where T\ = T2, it has been 
found to behave as a pow er-law in R with an exponent 
[i ~ —0.6 (Pan et al.ll2011l ). Unlike that case, the density 
is not a power-law in R when t\ ^12, instead having a 
well defined limit as R — > 0. Further, the clustering de- 
creases as e increases, as would be expected: the position 
of particles with significantly different stopping times are 
less correlated than those of particles with similar stop- 
ping times. 

From Figure [3] it is clear that the clustering diverges as 
e — > 0, as has been previously found for e = 0. However, 
because the clustering scales more weakly than eT 1 it is 
also clear, when combining that result with Figure^ that 
a net collision rate estimate um (f)/°( e ) does not diverge 
as e -> 0, which was also found in Paper I. Nonetheless, 



the confirmation of particle clustering even for different 
(but similar) stopping times continues to motivate the 
consideration of physics relevant to mostly-collisionless 
clusters of dust grains, which could otherwise have been 
dismissed as a result specific to the monodisperse case. 

It should be noted that the relatively weak power-law 
dependency on e for the clustering implies that the ex- 
treme levels of clustering predicted will not generate high 
dust-to-gas mass ratio clumps. For there to be significant 
direct back-reaction of the dust on the gas flow, the dust- 
to-gas mass ratio needs to be near unity. If the volume 
averaged dust-to-gas ratio is the canonical 0.01, this im- 
plies a 100-fold local increase in the dust density. To get 
a density enhancement over background of 100 through 
turbulent clustering, we estimate a maximum 6e range of 
0.29 <5e~ ' 77 = 100, or 5e — 5 x 10 -4 , or a maximum mass 
ratio ofl.5xl0~ 3 . As long as the dust is not tightly con- 
centrated at a specific grain size, the enhancement to the 
dust-to-gas mass ratio from turbulent clustering will be 
well below 100, and we do not have any reason to expect 
turbulent concentration to directly alter the character 
of the turbulence, meaning that our application of syn- 
thetic turbulence with no back-reaction possible is not 
immediately self- inconsistent. The streaming instability 
beco mes significant for dus t-to-gas mass ratios of about 
0.02 (|Johansen et al.l [200 91. double the canonical disk- 
averaged value. This would be expected for Se < 0.08, or 
a maximum mass ratio of 1.25. Triggering the streaming 
instability through turbulent concentration may there- 
fore be possible, but it will require almost all the dust 
being found in grains of very nearly the same size. 

4. VELOCITY DISTRIBUTION FIT 

It is now well known that the outcome of a collision be- 
tween tw o dust grains can de pend on the velocity of the 
collision (|Giittler et al.l [20Toh . This means that studies 
of collisional dust growth in protoplanetary disks need 
to include estimates of the rates and c ollisional velocity 
proba bilitie s, a step recen t ly be gun by iWindmark et al.1 
(|2012l ) and iGaraud et al.1 (|2012T ). However, those stud- 
ies assumed Maxwellian collisional velocity distributions, 
which is appropriate only for uncorrelated motion. The 
clustering data of Figure [3] immediately implies that for 
modest e, a Maxwellian fit is inappropriate. 

In Figure [4[ top panel, we show the normalized particle 
pair densities p(Ro/A, u) as a function of u p for our runs. 
While Rokp is not a constant, and so we are not cor- 
rectly normalizing our lengthscales to recover perfectly 
scale free behavior, we can nonetheless see a change in 
behavior as the stopping time ratio t\/t2 increases, go- 
ing from a sharply peaked distribution around a rela- 
tively low velocity to a much broader distribution around 
a significantly higher peak. We also have a qualitative 
change in the behavior of the high velocity tails, and the 
existence or non-existence of a low velocity bulge, which 
occurs between e = 0.5 and e = 0.9. Already by e > 2 
the high velocity tail appears to converge. 

As alluded to in Section 12.51 if a particle pair has rel- 
ative velocity u, then they can close a distance I ~ t\u. 
This means that the requirement of taking limits as 
R —> is less stringent for higher relative velocities. We 
illustrate this in the bottom panel of Figure [4[ where we 
show results from Run B for R — i?o/10 and i?o, along 




Fig. 4. — Top panel: The particle pair normalized densities 
p(Ro/4,u) as a function of u p . From lowest peak to highest we 
have runs A-K. There is a qualitative change as e grows above 1. 
Bottom panel: ~p(u) (black/solid), p(Rq/10,u) (red/dotted) and 
p(Ro,u) (blue/dashed) for run B. We can see that one needs to 
consider low values of R when evaluating p at low u, while larger 
values of R are both allowed, and required to sample an adequate 
number of particle-particle pairs when evaluating p for large u. 



with 



p(u) 



{p(Ro/10,u) u/u p < 0.5 

p(R /4,u) 0.5 < u/up < 1 

p(R /2,u) Ku/u p <1.5 

p(Ro, u) 1.5 < u/up. 



(7) 



This choice allows us to match the need for small R for 
low relative velocities while maximizing the particle pair 
count for rare, high velocity pairs and avoiding sharp 
discontinuities. In what follows we always refer to p. 

4.1. Limiting cases 

Before attempting to fit the entire distribution for all 
our stopping time ratios with a universal fit formula, we 
first consider subsections of the distribution. 

4.1.1. Low relative velocity limit 
In the limit of low relative velocity u, we find that 



p(u) 



for both large and small e, as shown in Fig- 



ure [S] This growth is faster for smaller e, which is reason- 
able as those particle pairs are more tightly correlated, 
and will collide at lower velocities. Note that the in- 
creased clustering for small e both allows us and requires 
us to consider smaller values of R than for large e. 

4.1.2. Large e - Maxwellian 

In the limit of large e > 1.25, the distribution of parti- 
cle pairs is well fit by an equation of the form 



p(u) ~ A 



(B + l)Cu 2 
B + u 2 




u [uj 




0.00 0.05 0.10 0.15 
u [u p ] 



0.0 0.016 0.033 0.05 
u [uj 



FlG. 5. — Low u limit of p(u) for four values of e. All are 
approximately quadratic in u. 




rz z 
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(8) 



Fig. 6. — Black/solid: p(u) for high e. Red/dashed: fits using 
Equation (01. 

as seen in Figure |51 where the fit function are the red 
dotted curves. However for e < 1.25 this shape no longer 
fits the results, as the tail flattens, approaching a lin- 
ear exponential. Equation is not quite a Maxwellian 
probability distribution, but it is close: for e > 1.25 we 
have B ~ 1, so the fit is nearly Maxwellian except for 
the high velocity tail. Note also that Equation has a 
strict maximum which means that our tail must eventu- 
ally hit zero. 

4.1.3. Low e, high u - Exponential 
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Fig. 7. — Black/solid: p(u) for low e. Red/dashed: fits consid- 
ering only the high velocity tail, using Equation {9}. 

In the limit of high relative velocity u, we find 

p{u) ~ ABe- Bu (9) 

is a good fit for low values of e as shown in Figure [3 
although our data range is too limited to constrain a 
polynomial prefactor. We show fits for e ranging from 
0.0625 to 0.9, and by the latter value the tail is clearly 
developing some curvature. 

4.1.4. Low e, low u - Clustering 

In Figure [5] we show low u zooms of p for our smallest 
e runs, along with exponential fits for the higher u tail. 
We can see that there is a distinct bump at low velocities 
that gets more pronounced as e gets smaller. In the e — > 
limit, this is the "cold" population of Paper I. Unlike that 
case however, for e ^ this bump in p is collisional as 
the peak collisional velocity um is located in the bump 
and is finite (see Figure [2] bottom panel). We attribute 
this enhancement to the particle pair number density at 
low u to turbulent clustering, and therefore refer to the 
population as the "clustered" population. 





4.2. Global fit - Form 




TABLE 2 
Fitting formulas 


Behavior 


Fitting formula 


quasi-Maxwellian 


fu = GM^f (26 MC 2 M + 1) j^e-W 


quasi-exponential 


/ B = a BCE (^# + l) j^re-— 


clustered 





IQ. 



rs; 10 
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Fig. 8. — Black/solid: p(u) for low e. Red/dashed: fits consider- 
ing only the high velocity tail, using Equation J9}. The difference 
between the two is evidence of the clustered population. Note that 
we have zoomed in to lower u, and the zoom is stronger for lower 
e due to the stronger, lower velocity clustering. 
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Fig. 9. — Normalized errors for our fitting formula ([f(u) — 
p(u)]/p(u)) using formulas from Table [2] The coefficients are de- 
termined for each run independently (listed in Table [3}- 

In this section we find a fit / for the form of p so 
that /(it) ~ p(u), that applies for all the e we consider. 
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Note that in this section, we consider only the functional 
form of /, and allow it to depend on coefficients that 
are determined for each value of e individually. We will 
fit the dependence of the coefficients on e in the next 
section. 

Combining the low e results, the high e results, and 
Figure SI we can see that there is a qualitative change in 
the behavior around e = 1, which corresponds to nearly a 
full order of magnitude in the mass ratio of the particles. 
This is interesting as by that scale particle-particle colli- 
sions can no longer be considered to occur between near- 
equals, and the larger grain has an increased chance of 
surviving even fairly energetic collisions. This change in 
behavior can be understood through Figure [2j a particle 
pair with relative velocity u samples a maximum separa- 
tion of order t\u. Particle pairs with e < 1 see low veloc- 
ity collisions that sample volumes well contained within 
a single eddy (or in our synthetic turbulence, a single 
gas velocity wavelength) , while higher e pairs can sample 
multiple eddies. This implies that analytical approaches 
to particle collisional velocities can only be applied to 
large e values: if particles begin their final approach when 
well separated, they do not start strongly correlated and 
an ensemble averaging approach can be used. If their fi- 
nal approach begins at a small separation however, then 
the initial conditions are critical and contain long time 
correlations which cannot yet be analytically treated. 

Our global fit model tries to capture the quasi- 
Maxwellian behavior at high e in a term /a/, the quasi- 
exponential high-w tails at low e in a term fg and the 
low velocity cluster for low e in a term fc, so our fit f(u) 
for ~p{u) is made up of 



f(u) = f M (u) + f E (u) + f c (u). 



(10) 



We will refer to the populations captured by the terms 
in order as quasi-Maxellian, quasi-exponential and clus- 
tered. The forms for Jm, Je and fc are shown in Tabled 
where the a parameters are amplitudes (note that the 
quasi-Maxwellian and quasi-exponential tail components 
are imperfectly normalized), the b parameters are offsets 
and the c parameters are inverse velocity scales. While 
the offset parameters have proven important to the fits, 
and both the quasi-Maxwellian and quasi-exponential be- 
havior is seen, eight fitting parameters allow for too much 
freedom: fitting the parameters for Runs A-K (using 
starting values from Tabled, we find values tabulated 
in Table [3l which shows some trends, but is not gener- 
ally directly useful (for example, the scatter in and 
&e is large). The fit quality is good however, as seen in 
Figure [HI where we plot [f(u) — p(u)] /p(u) . 

4.3. Master fit 

The values for the coefficients in Table |3] do not vary 
smoothly and monotonically with e, which implies that 
our fit formula is not a perfect match for the data. 
Nonetheless, we can construct a fit for the coefficients 
as a function of e which matches the data acceptably. 
Such a fit can then be adapted to generate an approx- 
imate p(u) for any e, and so can be used in a particle 
growth model. 

Requiring that for e < 0.5, cc = 2/um(c) which forces 
the peak of the clustered population to occur for u = Um, 
we can find fit formula in terms of e alone, arriving at 
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Fig. 10. — Normalized errors for our fitting formula ([/(w) — 
p(u)] /p(u)) using formulas from Table [2] with coefficients deter- 
mined by the formulas in Table [4] Note that e = 0.5 is shown 
with both the lower e fit (left) and upper fit (right). Compare 
with Figure [9] where the fit coefficients are determined for each 
run independently. 

Table |4] (certainly not unique), which identifies three 
regimes: a large e regime where the clustering term is 
insignificant, and two small e regimes where it is impor- 
tant. Significantly, a,E,bE and eg do not vary for small 
e: in this e range the large u tail is controlled by the ex- 
ponential terms, and they have converged as a function 
of stopping time ratio. 

We show the error of this fit in Figure QUI where we 
plot [f(u) — p(it)]/p(w), which is, unsurprisingly, inferior 
to that in Figure H] because the latter has full freedom 
for the coefficients of /. While for e = 1.25, 1.4, we 
underestimate cm , it is clear that for larger e our formula 
for cm does not approach the unity found in Table [3] 
rapidly enough. Less obviously, the formula for /m is 
imperfectly normalized. We know from Figure [3] that 
the integral of f(u) should be approximately 1 in the 
large e regime. For e < 2.5, we have J fM{u)du < <Xm, 
but for larger e it rapidly approaches its asymptotic value 
of cim ■ Those interested in such e values should consider 
normalizing /m directly. 

5. CONCLUSIONS 

We have extended the work of Paper I to find turbu- 
lence induced collision rates and velocities between parti- 
cles with different stopping times (i.e. of different sizes). 
Broadly, we find that the collisional velocity probability 
distribution is approximately Maxwellian for large stop- 
ping time ratios, but deviates strongly from that limit for 
particles with stopping time ratios less than 2 (which cor- 
responds to mass ratios below 10). We also find, similarly 
to Paper I, that our c ollisional velocities are smaller than 
previously estimated ([Ormel fc Cuzzill2007| ) for particles 
with modest stopping time ratios, although the previous 
estimates are not bad when the particles are significantly 
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TABLE 3 
Parameter fits 



e 




&M 


CM 


CLE 


b E 


CE 


«C 


cc 


0.0625 


0.53 


0.0013 


2.9 


0.35 


0.0016 


4.2 


2.2 


0.066 


0.125 


2.4 


0.0088 


5.2 


0.049 


13 


5.3 


0.59 


0.19 


0.25 


1.8 


0.051 


3.0 


0.11 


9.0 


5.6 


0.41 


0.21 


0.375 


1.8 


0.083 


2.4 


0.68 


0.12 


4.9 





N/A 


0.5 


1.2 


0.51 


2.2 


0.51 


4.0 


6.1 





N/A 


0.9 


1.6 


0.43 


1.3 


0.12 


6.1 


4.2 





N/A 


1.25 


1.6 


0.41 


1.2 





N/A 


N/A 





N/A 


1.4 


1.6 


0.71 


1.2 





N/A 


N/A 





N/A 


2 


1.5 


1.0 


1.0 





N/A 


N/A 





N/A 


2.5 


1.4 


1.7 


1.0 





N/A 


N/A 





N/A 


3 


1.2 


3.1 


1.0 





N/A 


N/A 





N/A 



TABLE 4 
Parameter function forms 



e range 


a M 


6m 


CM 


a E 


b E 




a-C 


cc 


e > 0.5 


1.6 


0.45e 


1 + 3^ 


0.04i d 


7 


6 





N/A 


0.25 < e < 0.5 


1.6 


0.3e 


6 

l+4e 


0.5 


0.33 


1.6 


0.08 
e 


2+1. 68e 

1.06e 


e < 0.25 


1.6 


4.8e 3 


6 

l+4e 


0.5 


0.33 


4.6 


0.08 

e 


2+1. 68t 

1.06e 



different in size. 

Further, we have managed to construct a usable esti- 
mate for the collision rates and velocities between dust 
grains of different sizes in turbulent flows which captures 
their qualitative behavior for arbitrary stopping time ra- 
tio (see Paper I for the case of monodisperse particles) . 

We interpret our results by noting that particle pairs 
with stopping time ratios close to unity (so particles of 
similar sizes) collide at low speeds. This means that they 
started their "final approach" while close together. The 
change in behavior between quasi-Maxwellian and non- 
Maxwellian collisional velocity probability distributions 
occurs approximately when particle pairs start their fi- 
nal approach one eddy size apart. This is the difference 
between particles whose positions and velocities have 
been correlated for a significant time, and particles which 
are nearly uncorrelated. Under that interpretation, it is 
hardly surprising that the collisional velocity probability 
distributions are different, nor that the large stopping 
time ratio collisions are quasi-Maxwellian. 

We find a counterpart to the cold, collisionless popula- 
tion of Paper I (collisions between particles of identical 
stopping times) , in our clustered population. In the case 
of non-unity stopping time ratios however, this popula- 
tion is collisional: the finite difference in stopping times 
forbids perfect correlation between the dust grains. We 
also find that while particles of similar stopping times are 



strongly concentrated by the turbulence, this clustering 
is unlikely to generate local dust grain density enhance- 
ments high enough above background to lead to stream- 
ing instabilities, although it is within the realm of the 
possible. 

While the different collisional velocity probability dis- 
tributions for different particle size ratios is not surpris- 
ing, it is of particular interest for dust coagulation mod- 
els. This is because the only effect that changing the 
velocity scale of collisions would have is to shift the final 
dust grain size, and shifting the final grain size enough 
that grains grow large enough to trigger gravitational in- 
stabilities or similar is implausible. Instead we want to 
enable more lucky dust grains to sneak through the frag- 
mentation barrier by a fortuitous sequence of low velocity 
collisions before growing large enough to sweep up small 
grains in high velocity encounters. 
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APPENDIX 

MEASURING DUST CONCENTRATIONS AND COLLISIONS 

In protoplanetary disks, dust grains that interact hydrodynamically with turbulence are far smaller than the s mallest 
turbulent eddies, unlike the sit uation in many other systems including rain droplets in the Earth's atmosphere (jShawl 
2003: Xu fc Bodenschatzll2008D . This makes it impossible to simultaneously resolve individual particle-particle colli- 
sions along with the turbulence that stirs them. This complicates the interpretation of numerical simulations designed 
to study the rates and velocities of turbulence induced particle-particle collisions. Rather than actually tracking in- 
dividual collisions, we must instead take periodic snap-shots of the entire system, identify interesting particle-particle 
pairs and convert their number, separation and relative velocities into global collision rates and collisional velocities. 
On the other hand, the same small size of the dust grains along with their high material density compared to the gas 
means that we can neglect the volume displaced by the grains, which is a significant simplification. 

To convert snap-shots of the system to collisional diagnostics, we must make the fundamental assumption that for 
any relative velocity it, there is a size scale id, simultaneously much smaller than the smallest turbulence length i v and 
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much larger than the dust grain radii a, where in all spheres of radius R < id constructed around test dust grains, the 
number density of dust grains with velocity u with respect to the test grains is independent of R. This assumption 
will be broken if new physics emerges when dust grains approach each other closely (such as mutual repulsion of 
negatively charged dust grains; lOkuzumil[2009h . If dust grains are large enough to induce wakes in the gas that affect 
particle-particle collisions, our assumption will also break down. 

Parameters and considerations 
The particles (moving with velocity u) and the gas (with velocity v) interact through the drag term 

where r p is the friction stopping time. Assume a turbulent cascade where the velocity of the turbulence at length-scale 
I is V£ and the turbulent correlation time (turnover time) at that scale is tt = l/vt . While it has not yet been rigorously 
determined, the strong expectation is that the turbulence which is most strongly coupled to pa rticle-particle rela tive 
motion (as opposed to bulk turbulent transport) is the turbulence at the scale £ p where tt = t p (jVolk et al.lll980| ). In 
the case of pairs of particles with different stopping times it is generally believed that the largest stopping time defines 
the most-coupled turbulence. 



Sampling method 

As in Paper I and this work, take snap-shots of the particles with temporal separation St, and bin particle pairs in 
space and relative velocity (velocity bin width Su). From the binned data calculate N(R,u,t), which is the number 
of particle pairs p\ and p 2 with \x% — x 2 \ < R, |«i — 1*2! € u ± 5it/2, in the snap-shot taken at time t. A dropped u 
implies summation over velocity bins while a dropped t implies averaging over snap-shots. 

Note that 



N(R,u,t)/ 







\ n v x 







(A2) 



is the effective number density of particles at relative velocity u ± Su/2 in spheres of radius R as seen by particles in 
the snapshot t. The volume averaged number density of particles po in a numerical simulation is however arbitrary 
(especially if they do not back-react on the gas), so it, along with the also arbitrary size of the velocity bin, should be 
scaled out: 



p{R,u,t) = N(R,u,t)/ 



1 



47T 



-i? J Po 



/Su 



(A3) 



(see also Equation ([5]), this paper) measures the density of targets with relative velocity u within distance R compared 
to the background number density. Similarly to N, for p a dropped u implies integration over velocity (giving the ratio 
of perceived number density to background) while a dropped t implies a time average over snapshots. A dropped R 
will be discussed in further detail below, but implies a limit R — > 0. 



Sampling rate 

If the relative velocity of a particle pair is u and 2R/St > u, then that pair can appear in multiple consecutive 
snap-shots. This can result in over/undersampling particles in regions of particle over/underdensities. If we have a 
large number of snap-shots taken over a total time interval t 3> R/u, our fundamental assumption implies that the 
oversampling of regions of over/underdensities will cancel. Nonetheless, it is safer to choose a sampling rate with 
St > 2R max /u m i n where R m ax is the largest particle-particle separation we consider, and v m in is the smallest relative 
velocity we consider. 

More importantly, turbulence generates significant structure in the particle distribution, which can take multiple 
particle stopping times or turbulent turn-overs to saturate, and leads to large swings in N(t), as seen in Paper I. As 
such, it is crucial to sample across multiple particle stopping times and multiple turbulent turnovers of the largest 
eddies which significantly interact with the dust grains. 



Sampling radii 

If R/u > T p , then the drag force will significantly affect the relative motion of particle pairs with a separation R 
before they could collide. This implies that for any collisional sphere radius R considered, there is a minimum velocity 
u > R/t p for which one can extract collisional information. Given that particle relative velocities are expected to 
couple most strongly to turbulence with turnover time tt p ~ t v , this immediately implies that a collisional sphere of 
the same size scale of that turbulence, l v = vt tt , will not be able to resolve collisional velocities of vt or less. This 
is a strong constraint: to calculate collision rates and collision velocities (especially modest velocities) for dust grains, 
we need to resolve particle separations much smaller than the length scale associated with the turbulence expected to 
maximally couple to the particle-particle relative motions. See the discussion of the "cold" population in Paper I for 
a particularly interesting consequence for particles with equal stopping times. 
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METHOD 

Taking the above considerations into account, one arrives at the procedure for extracting collisional diagnostics from 
the snap-shots. 

• Examine the time series of N(R,u,t) to verify that one is sampling a particle population whose collisional 
parameters have achieved a steady state, and that the sampling window is long enough to average over long-time 
fluctuations. 

• Take the limit p{u±) = limj?_> p(R, it*) for all values of it where this is possible, taking into account the require- 
ment that it* 3> Rmin/Tp. Where this is possible, the assumption is that one has achieved R < Id- Where it is 
not possible, either one has R > Id or the length scale £d does not exist. In the latter case, it is not clear how 
collisional diagnostics could be designed without resolving individual particle-particle collisions. 

Concentrations and rates 

• pop(R) is the number density of particles within spheres of radius R seen by test particles. If p(R) > 1 particles 
are clustered, while if p(R) < 1 they are segregated at a lengthscale R. 

• pop(u+) is the effective number density of particles with relative velocity u* that interact with test grains. Note 
both that it has units of density times inverse velocity and that it is not the term to use for collision rates! 

• o~pQU*p{u-k) is the contribution of the collisional velocity u* to the collision rate for a single particle, where a 
is the particle-particle collisional cross-section. Note the factor of it*, which accounts for the fact that particles 
moving slowly with respect to one another will collide only on long time scales. 

• apo J du u-kpiu-k) is the total collision rate for a single particle 

1/2 

• { \J du it*/?(M*)] / \J du it*p(ii*)l } is the root-mean-squared particle-particle collisional velocity. 
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